Fermions without fermi fields 



o 
o 

00 



X3 
O 

o 



> 
in 

oo : 
^ . 
a\ ■ 
o • 
^ : 
o . 

■ 

a: 

o 
o 



X 



R.C. Ball* 

Department of Physics, University of Warwick, Coventry CV4 7AL, UK 

It is shown that an arbitrary Fermion hopping hamiltonian can be mapped into a system with no 
fermion fields, generalising an earlier model by Levin & Wen[l, 2]. All operators in the hamiltonian of 
the resulting description commute (rather than anticommute) when acting at different sites, despite 
the system having excitations obeying Fermi statistics. Whilst extra conserved degrees of freedom 
are introduced, they are all locally identified in the representation obtained. The same methods 
apply to Majorana (half) fermions, which for cartesian lattices mitigate the Fermion Doubling 
Problem. The generality of these results suggests that the observation of Fermion excitations in 
nature does not demand that anticommuting Fermion fields be fundamental. 



As fundamental entities, fermion fields appear in con- 
fiict with the principle of locality: all fermion creation 
and annihilation operators anticommute no matter how 
far apart are the points in space at which they act 
(without restriction by causal connection). This fea- 
ture is built into Quantum Field Theory through the 
use of anticommuting Grassman fields, and supersym- 
metric string theories likewise have explicit anticommut- 
ing coordinates. Locality is only preserved in the Physics 
by conservation of fermion number, forcing the fermion 
operators to appear in pairs. In recent years there has 
been sustained interest in how particle statistics can be 
manipulated[3], and understanding the Fractional Quan- 
tum Hall Effect [4] has widenned the appreciation that the 
statistics of the elementary excitations of a system need 
not simply refiect the statistics of its components. Thus 
we ask whether it is really necessary to put fermion fields 
into physics 'by hand', or whether they can always be un- 
derstood as excitations emerging from quantum systems 
built of operators whose action is strictly local. 

It is well known how fermions relate to hard core 
bosons in one dimension[5] and how they can be built 
out of bosons in two dimensions with attached mag- 
netic fiux[6], and that neither approach extends natu- 
rally to three dimensions of space. The present letter 
builds on the recent work of Levin & Wen[l, 2], who 
showed that particular models of pairwise fermion hop- 
ping could be represented in terms of operators obey- 
ing locality, whilst the elementary excitations remained 
strictly fermionic. Their mechanism is essentially the re- 
verse of how Kitaev[7] showed that a particular hexag- 
onal lattice spin model has fermionic excitations. They 
showed that their mechanism worked on simple hypercu- 
bic lattices, in two and three dimensions explicitly, and 
interpret it in terms of string-net condensation[2, 8]. 

Here I show that fermion hopping on an arbitrary 
graph can be mapped into the excitations of a Hamilto- 
nian devoid of explicit fermionic operators, in which all 
operators acting at different sites commute. The dimen- 
sionality of space which the graph might approximate is 



irrelevant to the mapping, which is sensitive only to local 
coordination numbers (presumed finite). This strongly 
suggests that quite arbitrary fermion excitation spectra 
can be represented, and perhaps understood, as arising 
from the excitation of systems whose fundamental oper- 
ators obey locality in their commutation properties. 
We start from a generic fermion hopping hamiltonian, 



Hhop = ^ cfViCi + ^ 



Cj tijCj 



(1) 



where the fermion field Ci has standard anticommutation 
properties {cj,Cj} = {c^iCj'} = and {cf,Cj} = 
the Vi are simple 'on-site' potentials (relative to the 
fermion chemical potential) and the tij = t*^ are simple 
'inter-site' hopping matrix elements. The connectivity of 
the graph (or lattice) is encoded by which elements tij 
are non-zero, but below we will need to explicitly restrict 
the sum over links < ij >to those cases. We will also 
exploit the gauge invariance of the hamiltonian, that two 
models related by t^-^"* — tj^''e*''j'= are equivalent (through 
adjustment of phase of the Ck) provided the relative phase 
factors multiply to unity around all closed loops. 

We now introduce new operators Sij = —Sji modulat- 
ing the hoppings across each link, such that these opera- 
tors commute with each other and the original fermions, 
and with eigenvalues — —Sji 
is then generalised to 



= ±1. The hamiltonian 
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H gauge — clVjCj i SijcfuijCj 

i <ij> 

+ ^ ^ 9ij..z SijSj,..S,zSzi- (2) 
<ij..z> 

Here Uij — is^jtij so that the first two terms recover 
the hopping hamiltonian (1) for a particular set of out- 
comes s°- — — s°j of the operators Sij. We add extra cou- 
plings gij..z to (products round) 'Wilson loops' of the Sij 
operators; setting gij..z/ {sijSj ..s^z^^i) sufficiently nega- 
tive ensures that only combinations of the Sij eigenvalues 
which are gauge transformations of the original hopping 
hamiltonian contribute to the low energy states of the 
new hamiltonian. Because the operators Sij commute 
with each other and the hamiltonian, they could be sep- 
arately diagonalised to their eigenvalues Sy , which are 
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constants of the motion, and wc then recover the original 
tight binding hamiltonian (up to gauge symmetry) . Thus 
the system does still have its original fermion excitations. 

The key to obtaining a representation with locality is 
now to factorise each link operator into a pair of Majo- 
rana [half-] fermions, 

Sjk = i nijkmkj 

where the Majorana half-fermion operators are Hermi- 
tian (for simplicity) and have anticommutation proper- 
ties {Tnjk,mji ^k'} = '^Sjj'Skk' as well as anticommut'mg 
with all the original standard fermion operators a. We 
then associate each new half-fermion with the site of its 
first index, motivating us to rewrite the hopping terms 
in the hamiltonian as bfjUjbji, where 

b{j 7Tl{j C{ . 

Crucially the new operators bij commute, [6^ , bi/j'] = 
and [bfj,bi'j'] = 0, when their left (or site) indices 
are unequal, i ^ i'. The loop terms can also be ex- 
pressed in terms of operators conforming to locality in 
this way. Regrouping the factors in each loop gives us 
SijSjkSk.--S,zSzi = Bi^zjBj^ikBkj.-.Bz^.i where 

Bj^ik — iTTljiTTljk' 

We can finally eliminate all fermionic notation by writing 
= cfci and hence express the hamiltonian as 

Hgauge = ^ V^i^i + ^ b+^Uijbji 
i <ij> 

+ ^ 9ij..zBi^ziBj^i...Bz^.i- (3) 

<ij..z> 

From their definitions it is trivial to check that any pair 
drawn from all the operators rii , bij , 6^ , Bijk appearing 
in the hamiltonian (3) commute when their first indices 
are distinct. As a result we can factor the overall Hilbert 
space (of wavefunctions) on which they act into a prod- 
uct of single site Hilbert spaces, and the only non-trivial 
action of each operator is within the corresponding single 
site Hilbert space. 

We now focus on each site separately and note how 
the required commutation properties can be explicitly 
constructed. First note that each Dirac fermion can be 
split into a pair of hermitian Majorana half- fermions, 

2ck=mko+imk-i, (4) 

in terms of which Uk = (imfcoTifc-i + 1) /2. The com- 
plete set of operator properties then required on site k 
are now just the Majorana anticommutation relations 

{rrikm, rUkn} = ^Smn, -1 < TO, 71 < Zk, 

where z^. is the coordination number of (i.e. number of 
links to) site k. These relations are obeyed by standard 



(Euclidean) 4x4 Dirac matrices for 2 + Zk < 5 and by 
their 2''x2* generalisations for 2+Zk < 2s-|-l. It is crucial 
that we do not require to represent fermion anticommu- 
tation between local fermion operators on different sites, 
because all the terms in the hamiltonian contain an even 
number of fermion factors from each site. Writing jkii) 
for the A;'th Dirac matrix acting in the Hilbert space of 
the i'th site, we then have all the operators in the bosonic 
hamiltonian (3) expressed in terms of these: 

= ^7o(j)7-i(j) + ^, 

bjk = ^Ik'ij) (7o(i) + i7-i(i)) , 

Bj,ik =iji'{jhk'{j)- (5) 

Here on site j, 1 < k'{j, k) < zj denotes the local index 
associated with its link to site k. 

The end result is that all the operators appearing in 
the Hamiltonian arc expressed in terms of local matrix 
operators, which in turn are all equivalent to bilinear 
combinations of Dirac matrices: these might loosely be 
termed 'spins'. Where the original hopping hamiltonian 
had links with explicit fermion hopping, the derived gen- 
eralisation of the LW model[l] has coupling between the 
local spins. Most importantly, the spin operators for dif- 
ferent sites commute - yet by construction the system 
still has the original fermionic excitations. 

The loop products of operators Bijk can be further 
simplified in terms of products of new operators 

Pi] = lj'(u]){i)li'(3,i){j) (6) 

in direct index correspondence with the original op- 
erators Sij. These new operators are not equivalent 
to the Sij operators: in particular two operators Pij 
anticommute if they have one site index in common. 

Locality of all operators has been gained at the ex- 
pense of enlarging the Hilbert space. The original hop- 
ping hamiltonian (1) acted on a Hilbert space of dimen- 
sion 2^ where N is the number of sites, equivalent to 
N qbits, one for each site. To construct the generalised 
Wen Hamiltonian (3 with 5) we first added qbits equal 
to the number of links, A. However in the final local 
representation we carried fewer non-commutations, and 
multiplying the dimension of all the local Hilbert spaces 
leads to a dimension equivalent to total qbit count 

C = Ne + \no + a. 

Here Ne and Nq are the counts of sites with even and 
odd coordination number respectively, the even being less 
efficiently represented because the number of anticom- 
muting Dirac matrices is naturally odd. 

We should now expect that there are A — \No con- 
stants of the motion, and these and one more can be 
found as follows. First from every Hnk ij, we have the cor- 
responding Pij commuting with every term in the hamil- 
tonian (3). The result is that arbitrary (product) strings 
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of P operators commute provided they have no ends in 
common, including closed strings which have no ends. 
Secondly, for every even coordinated site we have her- 
mitian Tk = i*-^^^'^^ O^^-i Tj(^) anticommuting with 
every 7 on that site, and hence commuting with every 
term in the liamiltonian; however anticommutes with 
any P-string ending on site k. The maximal commuting 
set of all these operators then appears to be the union of: 

(a) all mutually inequivalent closed loop P-strings, equiv- 
alent to a minimal set of independent loop terms in the 
hamiltonian (3) , numbering A — N + 1 ; 

(b) the Ffe for all even sites, numbering iV^; 

(c) open P-strings ending on odd sites, with no ends 
in common and inequivalent under products with closed 
loop strings, all of which enumerate to by taking 
taking all the odd sites in (arbitrary) disjoint pairs. 

All of the above commute with each other and with the 
Hamiltonian, so we have in total A — \No + 1 explicit 
constants of the motion. Each corresponding operator 
has eigenvalues ±1, so its conservation removes one qbit 
from the dynamics. The author conjectures that the one 
extra conservation law relates to conservation of fermion 
number (modulo 2). 

Can one exploit the constants of the motion to reduce 
the size of the (quantum-mechanically coupled) Hilbert 
space? Trying this with Pij operators induces succes- 
sively less local anticommutation relations, tending to re- 
build the original fermionic representation. However we 
can ehminate the site- wise local Vi on even sites, leading 
us to a dimension-independent generalisation of what in 
one dimension corresponds to the Anisotropic Heisenberg 
Model representation of fermions. 

Let us focus on some particular even site i and in the 
following drop reference to that index. We can elim- 
inate 7-1 = i'-^"'"^''/^7o7i..72r, and then the operators 
in the Hamilonian which contained 7_i take the forms 
n = -4^7i-7.r + 1/2, 6fe = l7fe,7o (l - V-'^-n-.-fj:) , 
K = 57o7fc' (1 + i'^^^li-Tzr). Then because T still (evi- 
dently) commutes with the hamiltonian we can focus on 
the sector with eigenvalue F = 1 and make this replace- 
ment. Now the only matrices appearing are 7^, k = 0...Z, 
and we can take a minimal anticommuting representation 
of these[10], in terms of which we obtain n = ^ {1 + 70), 
bk = Ik'n, 6^ = n7fe' = 7/c' (1 — n) at each even coordi- 
nated site. 

For an arbitrary graph of even coordinated sites, we 
now have hamiltonian 

i <ij> 

+ Yl 9ij..zP.,P,.-P.zPzi, (7) 

<ii..z> 

and the explicit form for the Pij remains as per eqs. 
(6), except that the dimension of the local Dirac ma- 
trices has been halved. It is now particularly clear how 
the hopping terms conserve fermion numbers, and in- 
deed one can further show that nijj'{i)uijji'{j)nj = 



rii (1 — Uj) 7j/ (i)uijjii (j) (1 — Hi) fij making the inter- 
connection between occupation numbers IjOi and Ojl, 
totally expHcit. The loop terms are also particle con- 
serving and each P^ factor can be reorganised in similar 
manner. If one specialises to z = 2 corresponding to a 
one dimensional chain of sites, the Dirac matrices reduce 
to Pauli matrices and the first two terms of the hamilto- 
nian are exactly equivalent to the Anisotropic Heisenberg 
spin chain[9], well known to represent fermions in one di- 
mension. The periodic case has one loop term, but this 
reduces to a fixed scalar. It is gratifying that such a long- 
known fermion hamiltonian turns out to have natural 
extension to arbitrary dimensions and even to arbitrary 
graphs (of even coordination number). 

All of our analysis can be generaHsed to the case where 
we start from hopping of an arbitrary number h of half- 
fermions on each site. Hitherto we started from one stan- 
dard fermion per site in the fermion hopping hamiltonian 
(1) and split that into two half fermions (4), h = 2. In 
the general case any natural number /;, is allowed, with 
the onsite potentials and hopping matrix elements mak- 
ing arbitrary (hermitian) mixings amongst the h compo- 
nents. For h = 1 there are no on-site potential terms, 
and hermiticity restricts the hopping matrix elements to 
be pure real, Uij = u[j, giving us 

Hhop,i = i ^ u-j-moiTOoj. (8) 
<ii> 

The standard fermion case discussed eaher is /i = 2, and 
it reduces to the sum of two (commuting) h = 1 hamil- 
tonians if the onsite potential terms are zero and the Uij 
are all pure real (to within a gauge transformation) . The 
case /i = 3 commands interest for particle physics. 

The discussion of conservation laws and reduction of 
Hilbert space generalises quite trivially, with the under- 
standing that even and odd sites are identified by whether 
h + zis even or odd. Particularly simple spin models are 
obtained from h = 1, where the local operators bij are 
hermitian and can be directly represented as single lo- 
cal Dirac matrices 7j(?) (rather than bilinear products) 
requiring only Zi Dirac matrices at each site, and with 
Bijk i^j{i)^k{i)- Quite general Majorana fermion 
hopping hamiltonians can therefore emerge from the ex- 
citations of models built out of simple local spin opera- 
tors. The gauge symmetry changes naturally with h. For 
h = 1 we have only Z2 (sign changes), corresponding to 
the transformations exploited in the 'KLW trick' of intro- 
ducing the S operators, and for /i = 2 we already noted 
U{1) = 0{2) gauge symmetry. For /i = 3 we would have 
gauge symmetry 0(3) if the onsite terms are sufficiently 
symmetric, but quite arbitrary mixing of components al- 
lowed in the onsite terms of the hamiltonian would in 
general reduce the gauge symmetry to 0(2). 

The simple (hyper-)cubic lattices are of special interest 
as a source of insight into the continuum limit, and parti- 
cle physics interest focuses particularly on the case where 
the fermion excitations are (almost) massless. The spe- 
cial case of the main analysis above previously presented 



4 



by Levin and Wen, which was for simple square and cu- 
bic lattices with h = 2, falls into this category. However 
that work also suffered from the standard 'Fermion Dou- 
bling Problem' of lattice fermion models[ll], yielding 2'' 
massless fermions in d dimensions. 

Starting from a single half-fermion {h = 1) hopping 
model halves the Fermion Doubling to give just 2^*"^ 
massless fermions in d dimensions. From the point of 
view of the spin models, this is just as natural a starting 
point as h = 2. For simplicity we present the calcula- 
tion in one dimension, which suffices to demonstrate the 
novelty, as the elaboration to higher dimensions simply 
follows in the same manner as Levin & Wen[l, 2]. A 
simplest one dimensional Majorana hopping hamiltonian 
is 

Hi = -w'^m„m„+i 

n 

and imposing periodic boundary conditions for such a 
system of N sites results in a spectrum of standard mass- 
less fermions. 

Hi = 2u'Y^ {^ctck - l) sinfc. 

Here the wavevectors span Jialf the first Brillouin zone 
of the lattice, k — n2'K /N , n = 1 to N/2, the fermion 
operators are given by Ck = l/\/2]V5]]„ e~*'^"mn and 
the negative wavcvcctor Fourier components give their 
conjugates. These Ck are standard full fermion operators, 
in particular obeying {c^, c;} = Ski- 

The above is in effect a staggered fermion[12] solution 
to the Fermion Doubling Problem[ll]. The full fermions 
obtained by halving the Brillouin zone can be mapped 
onto one full fermion per two original sites, but these 
are non- locally related to the original half fermions. For 
example associating the full fermions with even sites gives 




2q + l 



This parallels another solution to the Fermion Doubling 
Problem, in which highly non-local hopping with am- 
plitudes similar to the above is used to force a better 
spectral approximation to the continuum[13]. Here we 
have local hopping for the Majorana particles and the 
non-local amplitudes encoded in how they map onto con- 
served full fermions. We are not forced to work with 
these non-local realtionships explicitly, even in the tight 
binding hamiltonian, because we can always make full 
fermions more locally, for example out of adjacent pairs 
of half fermions; the price of doing this is that in terms 
of the latter the hamiltonian then contains pair creation 
and annihilation terms. 

One can introduce mass without doubling the spec- 
trum, simpl y by m odulating the bond strengths 
V u'^ + + e(— 1)" to open up an energy 
gap at i? = (corresponding to A; = 0,7r). One 



then finds dispersion relation£'(/;;) = ±4 \/ 7/2 hIh^ + 
with two states at each wavevector over the quar- 
ter Brilliouin zone < k < 7r/2. Adding more 
cartesian dimensions with suitably frustrated signs to 
the couplings leads (as per Levin and Wen[l, 2] for 
the full fermion case) to the natural generalisations 
in higher dimensions, for example E{kx,ky,kz) = 

±4^^ u'^ sin^ kx + u'^ sin^ ky -\- uf sir? kz + in three di- 
mensions. Each increase in dimension generates a dou- 
bling of the spectrum, to two species or a two component 
fermion in d = 2, and a four component fermion in d = 3, 
matching standard free fermion fields. 

In overall summary, we can now construct hamilto- 
nians in terms of strictly local operators (in the sense 
of their commutation relations) whose excitations corre- 
spond to fermion hopping on whatever graph is desired. 
If one accepts inter-site hopping as an approximation to 
continuum particle dynamics, then this means that al- 
most any fermion problem can be constructed out of the 
excitations of what might loosely have been described as 
a Bose system. 

The method has been presented for free (non- 
interacting) fermion systems. However it is trivially gen- 
eralisable to the case of an arbitrary interaction, provided 
this is a function of the state number operators (or other 
bilinear combinations of local fermion variables) . In par- 
ticular Coulomb interactions between different sites are 
allowed, as are Hubbard interactions. Modifying hopping 
according to occupation numbers, as in t — J model, is 
also readily incorporated. 

The complete and relatively local enumeration of all 
the constants of the motion may have consequences for 
string-net condensate interpretations[2, 8]. By introduc- 
ing terms in the Hamiltonian coupled to all the conserved 
quantities, we are free to remove all the degeneracy of 
even the gauge-equivalent combinations of S link opera- 
tors introduced in eq. (2) from low energy states, in as 
many ways as the degeneracy introduced. If string-net 
condensate states survive this splitting, they have to be 
correspondingly degenerate to start with. 

The 'spin' hamiltonians constructed to give exactly 
the specified fermions are somewhat cumbersome. Are 
there simpler local spin hamiltonians which still give 
fermion excitations? The generalised Anisotropic Heisen- 
berg hamiltonian is one case in point: setting the onsite 
terms to zero and removing the number operator factors 
from the hopping terms turns out to give the same as 
starting from a Majorana hopping hamiltonian (8). Also, 
we could by deliberate construction build hamiltonians 
which are supersymmetric in their excitation spectrum, 
just by adding in the desired bosons. Are there variations 
on the theme which are more naturally supersymmetric? 
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